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Introduction. In the mathematical finite element modeling, an 
average value of the mechanical characteristics of the deform- 
able solid material is used. In aircraft, machine building, con- 
struction engineering, medicine and other fields, polymer 
composite materials and materials of natural origin are in- 
creasingly used. In the latter case, the actual change in the 
mechanical characteristics differs significantly from the aver- 
aged change; therefore, when using the averaged parameters to 
build and analyze finite element models, the results can be 
significantly distorted. This paper describes the creation of 
mathematical methods for studying changes in the mechanical 
characteristics of a material of inhomogeneous deformable 
solids. The results obtained in this way are used to construct 
finite element models and analyze their stress-strain state. 
Materials and Methods. Naturally occurring materials and 
composites are considered as inhomogeneous deformable 
solids. To study the changes in the mechanical characteristics 
of the material, a method was developed based on the use of 
two components: the pixel characteristics of raster images 
scanned by a computer tomograph and the experimental data 
of field tests of standard samples. 

Research Results. A complex of mathematical methods has 
been developed for modeling the interpretation of scanning 
raster images by a computer tomograph, which allows for the 
study of any complicated structures of real deformable solids. 
The results are used in the construction of finite element mod- 
els of such bodies considering the heterogeneity of the me- 
chanical characteristics of the material. The analysis of the 


stress-strain state of finite element models of test samples has 
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Beedenue. IIpu MaTeMaTHu4ecKOM KOHCYHOIJICEMCHTHOM MOJIe- 
JIMpOBaHHN UCIIOIIb3yeTCA ycpeqHeHHoOe 3Ha4eHHe MexaHHye- 
CKHUX XapaKTepHCTHK MaTepHasia WePOpMupyeMbIxX TBep/IbIX 
Tel. B aBiva-, MallIMHOCTpOeHHH, CTPONTeJIBCTBe, MeJIMIIMHe U 
Mpyrux oOsacTax Bce WIMpe NpMMeHAIOTCA MOJIMMepHbIe KOM- 
TO3HTHbIC MaTepHasibl UW MaTepHasIbI IPMPOAHOTO MpOUCxOx- 
yeHusa. B nocieqHeM ciryyae peasIbHoe H3MeHeHve MexaHnye- 
CKUX XapaKTepHCTHK 3HA4UMTeIBHO OTIMY4aeTCA OT ycpeyHeH- 
HOrO, CJIeqOBaTeJIbHO, pH UCHONb30OBaHHH yCpeHeHHBbIX 
TapaMeTpoB JIA MOCTpOeHHA UH aHasIM3a KOHCUHOSJIEMCHTHBIX 
MOJeJIeH pe3yIbTaTbI MOTYT CYI[eCTBeEHHO MCKaxKaTBCA. 

B WaHHOM cTaTbe OMMCaHO co3qaHWe MaTeMaTH4eCKHX MeTO- 
JOB UCCUeOBAaHHA W3MeCHCHUA MeXaHM4eCKUX XapaKTepHCTHK 
MaTepHasia HeEOJHOPOHbIX Je~OpMUpyeMbIxX TBePbIX Ted. 
IlostyueHHbIe TakKHM 0O0pa30M pe3yJIbTaTbI NPHMMeHeHbI WIA 
TOCTPOeCHHA KOHCUHOSJICEMCHTHBIX MOJIeIeH WM aHalIM3a UX 
Hallps9KeCHHO-JJe~OpMHPOBaHHOrO COCTOAHHA. 

Mamepuazoi u memooovi. B KayecTBe HeEOAHOPOAHBIX Wedop- 
MUPYyeMBIX TBeEPbIX TeJI PaCCMOTPpeHbI MaTepHasIbl MpUpor- 
HOrO IIpoOwcxoxKTeHHA WU KOMMO3HTBI. Jia UccueqOBaHHA 13- 
MeHeHHA MeXaHW4ecKHX XapaKTepHCTHK MaTepHasia pa3pavo- 
TaH CHOCcOO, OCHOBaHHBIM Ha UCIMONb30BaHHH JIBYX COCTAaBJLA- 
FOWIMX: MUKCeJIbHONW XapaKTepHCTHKH pacTpOBbIX H300paxe- 
HUH CKaHHpOBaHHA KOMIIbIOTePHbIM TOMOTpadomM HU 3SKCIIepu- 
MCHT@JIbHBIX JI@HHbIX HATYPHbIX UCMbITAaHHN CTAaHAapTHbIx 
oOOpa3I{OB. 

Pe3yibmamol ucciedoeanus. Co3aH KOMIMIeKC MaTeMaTH4ue- 
CKUX MeTOXOB MOJeMpOoBaHHA WHTeplpeTaluu pacTpOBbIx 
W300paxKeHHH CKaHHMpOBaHHA KOMIIbIOTepHbIM TOMOrpadom, 
NO3BOJAFOIIMH MpOBOAMTb UCceOBAaHve JOObIX CJIOXKHBIX 
CTPyKTyp peasIbHbIx JJePOpMUpyeMBIX TBepJIX Tes. Pesyb- 
TaTbI HCIOJb3yYIOTCA pH MWOCTPOeHHH KOHCUHOSIIEMCHTHBIX 
MOjleIeH TaKHX TeJI C yYeTOM HeEOJHOPOHOCTH MexaHnye- 
CKUX XapaKTepHCTHK MaTepHasia. 

AHasIv3 Hallps9KeHHO-eMOpMHpOBaHHOrO COCTOAHHA KOHeY- 
HOSJICEMCHTHBIX MOJleIeH TeCTOBBIX OOpa3IlOB OKa3ajI TOU- 
HOCTb HW CXOJJMMOCTS YMCIICHHOFO pellleHHA MeTOJa KOHCUHBIX 
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proved the accuracy and convergence of the numerical solu- 
tion of the finite element method in modeling the property of 
heterogeneity of the mechanical characteristics of the material. 
Discussion and Conclusions. The developed approach can be 
applied to any physical principles of scanning (X-ray, ultra- 
sound, laser, etc.) and for any types of materials if the data 
obtained as a result of scanning is developed in the form of a 


digital (raster) image. 


Keywords: finite element method, deformable solid, inhomo- 
geneity, mechanical properties of material. 
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9JIEMCHTOB TPH MOJeIMpOBaHHH CBOMCTBA HeEOJHOPOHOCTU 
MeXaHHYecKHXx XapaKTepHCTHK MaTepnasia. 

Oécyacdenue u 3akiouenuaA. Pa3spadoTaHHbI MOAXO, MOKET 
ObITh MPUMeHeH JIA WHOOBIX (U3NH4eCKHX MPHHUMMOB CKaHH- 
poBaHua (peHTTeHOBCKHH, yIbTpa3ByKOBOH, Ja3epHbI U Jp.) 
WM IA WKOObIX THMOB MaTepuaoB, ecu UHDOpMallnsA, Wosy- 
YeHHad B pe3yibTaTe CKaHUpOBaHHA, CPOPMUPOBAaHA B BUC 


uudpoBoro (pacTpoBoro) v300paxKeHHaA. 


Ku110ueBble CJIOBa: MeTO], KOHCYHBIX 3JICEMCHTOB, edOpMu- 
pyeMoe TBepoe TesIO, HEOJHOPOHOCTb, MexaHWyeckne xa- 
pakTepHCTHKU MaTepnalia. 
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MaTH4eCKOFO KOHCYHOSJICMEHTHOTO MOJIeIMpOBaHHA HeEOHO- 
POHBIX JE@OpMUpPyeMBbIX TBepJIbIX TeJI C IPHMeHeHHeM CKa- 
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Introduction. In the finite element (FE) modeling, mechanical characteristics of the material of deformable 
solids (DS) are specified as an averaged value used for the whole model. For example, mechanical characteristics of 
steel are determined through tensile and compression tests of standard regular samples. This approach is acceptable for 
the design and engineering practice. However, in modern aircraft, engineering, construction, medicine and other fields, 
polymer composites and materials of natural origin are increasingly used. A real change in the mechanical characteris- 
tics of such materials differs drastically from an average value; therefore, when using average parameters to build and 
analyze finite element models, the results can be significantly distorted. 

Research Results 

1. The main mathematical dependence of modeling heterogeneity of real deformable solids. To solve the 
presented problem, a method for determining a real change in the mechanical characteristics of the material based on 
the use of: 

- DS scan results; 

- data of field tests of standard samples [1 ]. 

This method makes it possible to recognize any complex heterogeneous material structure and apply this data 
to increase accuracy and realism in mathematical finite element modeling of real DS. 

Fig. la shows an example of a raster image of scanning a material of natural origin (human bone). Its section 
demonstrates a change in mechanical characteristics depending on the values of the color indices (pixels) [2]. Fig. 15 
shows the scan result of the composite material in the form of a raster image of the cross section of a regular sample 
(the scanning was performed by V.G. Tolstikov). Real changes in the mechanical characteristics of the fibers and matrix 
(adhesive) in the structure are indicated. In addition, the scanning has shown the possibility to assess quality of the 
composite materials (fiber arrangement, adhesive layer thickness, internal defects, etc.). 
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Fig. 1. Material of natural origin (a) and composite material (6) scan bitmap analysis 


A general diagram of mathematical simulation and the algorithm of the proposed scanning technology are 


shown in Fig. 2 [2]. 
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Fig. 2. Mathematical and finite element modeling of real deformable solids 


The diagram reflects a full range of modules: from scanning DS to obtaining its FE model and analysis of the 
stress-strain state (SSS). If necessary, the cycle of operations is ensured within each block of the algorithm. 

Scanning is performed using a computer tomograph (CT). When working with deformable objects, it provides 
image packetizing; each one is a raster image in a certain section [3, 4]. Standard raster image formats are jpeg, dicom, 
bmp, png, etc. Regardless of the format used, a raster image should be represented as a numerical matrix for further 
processing. 

To determine the dependence of mechanical characteristics (for example, in the form of an elastic modulus) on 
the values of color indices, a special weight factor is applied [1], which establishes the relationship between the average 
value of the color indices n,, and the experimental value Eg. The average value of the elastic modulus (£97) obtained 
in the experiment on the tension (compression) of standard samples is used to interpret mechanical changes in the pixel 
characteristic. 

Then the weight factor that determines the transition from the pixel color index to the elastic modulus is deter- 
mined through an expression of the following form: 

kg = 208 | (1) 
Nop 
Here, n,, is an average value of color indices of a scanning bitmap packet: 
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I, j) 
i=l j=l 
1.» = ——————_. (2) 
= mn 
However, for a large amount of CT SSS scan data with a high degree of the structure heterogeneity, it is re- 
quired to use the mathematical expectation dependence when determining an average value of the color index. That is, 


m on 
=! 


Np 18 defined as an average value: 
Ncp = E(n) = yo mPi (n;), (3) 
i=l 
where n; is the color index value on the interval J; p;(n;) is the probability of occurrence of n; (color index values) in the 
interval J. 

Pixels, each of which has a corresponding color index, are assigned the values of the elastic modulus of the 
material E;. For this, the weight factor kg is used: 

E, =n; kp. (4) 
Here, £; is the elastic modulus value corresponding to the value of the color index number n,; kz is weight factor of the 
elastic modulus. 

The dependence (4) for the FE model is implemented through solving the following tasks: 

- determination of the dependence of mechanical characteristics on the values of color indices; 

- determination of the mechanical characteristics of the material in the nodes and in the finite elements of the 
FE model. 

A change in the elastic modulus is determined on the basis of a linear or non-linear dependence on the values 
of the color indices. For a material with a low degree of heterogeneity, the dependence can be constructed according to 
the linear law (4) using a single constant value of the weight factor of the elastic modulus. 

However, such materials are rare in nature and engineering. Most often, several areas can be distinguished in 
the structure of a DS material, each of which has a different variation range of the elasticity modulus than the others. 

It is proposed to apply a nonlinear law for such materials, which provides its own weight factor value in each 
region of the DS structure. In this case, the dependence of the change in the DS elastic modulus on the color indices is 
presented as a piecewise linear function (Fig. 3a). 
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Fig. 3. Piecewise linear (a) and continuous (5) function of DS elasticity modulus 
variation from color index values 


Using a piecewise linear function is difficult from the point of view of the computational algorithm efficiency, 
in particular, due to the following reasons: 

- the larger the number of the DS structure regions, the greater the number of weight factors and the more 
complicated the calculation of the elastic modulus for nodes and finite elements; 

- transition values of color indices between DS regions or objects involve not one value, but ranges of values of 
color indices in a bitmap (see Fig. 35). 

To solve this problem, a nonlinear dependence in the form of a spline, for example, a cubic one [5], can be 
constructed on the basis of a piecewise linear function. Such a nonlinear dependence of the DS elastic modulus on color 
indices is universal when combining several zones with different patterns of mechanical characteristics (see Fig. 35). 
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Next, the interpretation of the scanning data is associated with the interpolation and transmission of color index 
values from pixels to the nodes of the FE mesh. In the most complex and general case, the FE nodes of the mesh lie 
between adjacent planes with a bitmap. They vary widely in shape and size. Besides, additional internal contours and 
other geometric features appear in them. In this case, a number of mathematical dependences are used for interpolation; 
they are determined with respect to a straight line inclined to the Z axis [6]. 

In [2], a more similar description of the mathematical dependence and the block diagram of the algorithm for 
determining the mechanical characteristics of the material in the nodes and finite elements of the DS FE model ts pre- 
sented (see Fig. 2, block 7). 

To use the results of determining the heterogeneity of the material mechanical characteristics in the FE model, 
they are converted to a file with the *pc/ extension. Two functions of the Patran Command Language: “material.create” 
and “elementprops_ create” [7, 8] are used for this. 

2. Investigation of the accuracy and convergence of the analysis results of the FE models considering 
simulation of heterogeneity of mechanical characteristics. In the framework of this study, accuracy and convergence 
of the numerical solution by the finite element method (FEM) using the presented technique are determined. For this 
purpose, we used the SSS calculations of a human bone tissue (a fragment of the femur). 

The selection of DS and its analysis are nonessential, but they are caused by important factors: heterogeneity of 
the bone material, individuality of its geometry, as well as high technology and quality of CT scanning in medicine. 
Bone tissue is well studied in practice, which makes it possible to reasonably consider the correctness of its mathemati- 
cal simulation. 

The first study uses parallelepiped specimens cut from a fragment of a human femur (Fig. 4a). The results of 
the study are shown in Fig. 45 as proof of accuracy and convergence of the FEM numerical solution with respect to the 
full-scale experiment data [9, 10]. 
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Fig. 4. Dimensions and orientation of real samples (a); graph of convergence of normal stress in the center of samples (5), 
longitudinally (1), circumferentially (2) and radially (3) 
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The results show that to achieve the desired accuracy of the FEM numerical solution, a FE mesh with a density 
of three or more finite elements per 1 mm is required [11]. It is also obvious that the heterogeneity property of the DS 
material mechanical characteristics can be reflected by a set of finite elements with their own elastic moduli and the 
isotropic structure of the material. 





a) 
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Fig. 5. SSS analysis result using HEX (a) and TET (bd) type finite elements; graph of convergence results (c) 


To study the applicability of two different types of solid FE (hexahedron and tetrahedron), an additional study 
was conducted using incompatible shape functions (see Fig. 5) [12]. The results show that hexahedron-type FE has the 
best indicators in terms of accuracy and convergence of results, as well as in the resource costs. 

3. Modeling of real deformable solids. The technology presented is applied to real objects. For this purpose, 
the SSS of the FE of tooth setups with a defect and a sealant made of composite material was constructed and analyzed. 
One model presents a defect in the form of a carious region (cavity) in the upper part of the enamel. In the second 
version of the FE model, this area is filled with a sealant made of composite material. Both models are built considering 
heterogeneity of the material mechanical characteristics and an individual geometry. 

A wireframe model of a human tooth, shown in Fig. 6a, is based on the scan results [13, 14]. A change in the 
material mechanical characteristics in the solid FE model (Fig. 65) is simulated in a special program [13, 14]. Input 
data: results of scanning and field testing. 
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Fig. 6. Construction and analysis of FE models of a human tooth considering heterogeneity of material mechanical characteristics: 
a) wireframe model; 5) change in elasticity tooth modulus in FE solid model; c) design scheme; d) solid geometric model; 
e) FE mesh generation; f) SSS analysis result of FE model of tooth with simulated defect; g) SSS analysis result of FE model of tooth 
with simulated composite fillings 


The design scheme is presented in Fig. 6c. A defect and a composite filling are modeled on the upper tooth 
(see Fig. 6d and 6e). 

The analysis result of the FE model with a defect (see Fig. 6f) shows that the carious part of the tooth does not 
absorb the load. That is, a healthy part of the tooth has an additional load, and the probability of decay increases. The 
analysis result of the second FE model, after filling the caries area with a sealant made of composite material (see Fig. 
6g), shows that the tooth almost completely restores its functional capacities since the stress pattern is identical to the 
picture obtained for a healthy tooth [2]. 

Fig. 6g shows that stress maxima are observed at the contact points on the tops of the dental tubercles. A 
higher load is picked up by the area of enamel and dentin. The inner part of the tooth is less stressed. This result is 
compared to the calculation [2]. An average value of the mechanical characteristics of enamel and dentin is used. 
Heterogeneity of the structure of the tooth material is not taken into account. The comparison shows a more uniform 
stress distribution in the tooth body when modeling the distribution heterogeneity of its mechanical characteristics. 


Discussion and Conclusions. The study results validate the accuracy and reliability of the presented method of 
FE modeling based on scanning real DS considering heterogeneity of the material mechanical characteristics. 
The presented complex of mathematical methods for modeling the interpretation of CT scanning bitmaps al- 


lows the study of any complex structures of real DS. The results of such simulation are used to construct FE models 
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considering heterogeneity of the mechanical characteristics of the material. 
The developed mathematical modeling technique can be applied for any physical scanning principles (X-ray, 


ultrasound, laser, etc.) and for any types of materials if the data received is a digital (raster) image. 
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